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1 Introduction 

Metastable phase decay on the several types of heterogeneous centers remains 
a rather actual problem for theoretical investigation. For the first time the 
theory for the kinetics description was constructed in p| . This approach de- 
composes the general situation into characteristic situations which are rather 
simple. All limit situations characterized by small values of characteristic 
parameters can be solved by the slightly modified versions of the iteration 
method initially proposed in Only the intermediate situation requires a 
special method of consideration which is based on the special monodisperse 
approximation. 

When the total number of the heterogeneous centers is one and the same 
for different types of centers then as it is shown in there are only two 
characteristic situations: the intermediate situation and the situation of the 
strong unsymmetry. 

The special monodisperse approximation is well based and can be spread 
on a more general situation. This generalization is important not only to 
present the description in a more compact form. The special monodisperse 
approximation allows to reduce the error appeared in the limit situations. 
Really, in the situation of the strong unsymmetry (this is a standard limit 
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situation) one has to use the monodisperse approximation. If we use the 
special monodisperse approximation instead of the already used one we shall 
reduce the error. 

The universal character of obtained solution was shown in 0]. 

The general recipe to use the special monodisperse approximation was 
suggested in ^ in the abstract manner. So, now it rather natural to show 
how the special monodisperse approximation works in the limit situations 
concretely. 

This publication is intended to show how to use the special monodisperse 
approximation in the limit situations. We shall use for example the situation 
of the strong unsymmetry which appeared when the number of the first type 
centers equals to the number of the second type centers (see [Q). This sit- 
uation together with the intermediate situation completes the consideration 
of the general case . 

2 Formulation of the problem 

Consider the system with two sorts of heterogeneous centers (they are marked 
by subscripts 1 and 2). Suppose that the total numbers of heterogeneous 
centers rjtoti and rjtot2 are equal 

vtoti = 'ntot2 (1) 

At the initial moment of time which is denoted by the subscript * there 
exists only heterogeneous centers. We shall call the period of intensive for- 
mation of the droplets on the heterogeneous centers of given type as the 
nucleation on the centers of given type. 

We shall use in the estimates some characteristic values. Denote by Ait 
the duration of the nucleation on the first type centers and by the dura- 
tion of the nucleation on the second type centers. 

Consider the situation when the rate of nucleation on the first type of 
heterogeneous centers strongly exceeds the rate of nucleation on the second 
type of heterogeneous centers. We shall choose the sorts of heterogeneous 
centers to have 

where /*i is the amplitude value of the stationary distribution function. 
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Due to (|1|) the last inequality is practically equivalent to 

exp(-AFi) > exp(-AF2) 

where AF is the height of the activation barrier. 

The power of met ast ability is characterized by the supersaturation 

( = (n - noo)/noo 

where n is the molecules number density of the vapor, Uoo is the molecule 
number density of the saturated vapor. 

Then the balance of the substance in the closed substance requires 

C = C + Gi + G2 

where Gi is the number of the molecules in the liquid phase around the first 
type centers (taken in units of rioo), G2 is the number of molecules in the 
droplets formed on the second type centers. 

We shall describe the droplet by the value of dimensionless radius p which 
equals to the cube root of the number u of the molecules inside the droplet 

p = u'/' 

This value is convenient because the rate of its growth is one and the same 
for all sizes 

dp 

where the constant r is the characteristic time. The last equation is valid for 
the supercritical droplets under the free molecular regime of the substance 
exchange. 

Then automatically one can see that the distribution /j(p, t) of the droplets 
over p in the moment t depends only on the intensity of the droplets forma- 
tion in the time when the droplet of given size were formed. 

The application of the ordinary monodisperse approximation ("total monodis- 
perse approximation" ) is based on the following explanation: 

- Suppose we can suggest an approximation 



Gi ~ NuotPl/noo 
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where Nuot is the total number of the droplets formed on the first type centers 
during the condensation process. 

The last approximation is good when po is many times grater than the 
value of Po at the end of the period of intensive formation of the droplets on 
the first type centers, i.e. po(Ai)f:): 

Po(t) > Po(Ait) 

This approximation is important when Gi stops the formation of the 
droplets on the second type centers, i.e. at /S.2t. So, it is necessary to be 

Po(Ait) < po(A2t) 

As far as dp/dt is rather smooth function of time one can rewrite the last 
estimate as 

Alt < Aat (2) 

One can estimate the number Nuot of the droplets formed on the centers 
of given sort as 

Nitot ~ Ji*^it 

where 

^2* — fi^^oo^/ C 

is the initial rate of nucleationQ. Then the violation of (0) means that 

Nltot > N2tot 

As the result the number of the droplets formed on the second type centers 
is negligible. 

The negligible value of A^2tot in the unique situation when the ordinary 
monodisperse approximation fails was the ground to apply this approxima- 
tion in the situation of the strong unsymmetry p[. But we see that if we are 
interested in the value of A^2tot without any respect to Nuot then the question 
is still open. Below we shall resolve this problem. 



^/i* is expressed in units of ric 
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3 Special monodisperse approximation 

For Gi one can write the following relation 

roo 

Gi= / p'Mp,t)dp 
Jo 

Now one can analyze the subintegral function 

9i{p,t) = P^Mp,t) 

which has the sense only for positive p. 
One can see that 

9i{p,t) = 

for 

p > po{t) ^ f ^dt' 

Jo T 

One can see that 

9i (P) ^) ^ P fi* — 9i appr 

Consider now the pseudo homogeneous situation. It means that we ne- 
glect the exhaustion of the heterogeneous centers by the droplets. Here 

9i{p,t) ~ P^fi* 

for all p from 

Po-p< (0.7-0.8)CA,Vt 

As far as p^fi* is the sharp function of p we see that gi is even more 
sharp function of p. Then it is quite reasonable to suggest the monodisperse 
approximation for g^. 

To construct the monodisperse approximation for gi one has to solve how 
to cut the tail at small p. Despite the rapid decrease at small p the tail can 
not be integrated at least on the base of approximation gi appr (if we forget 
about the restriction p > 0). There are two ways to do it. 

The first way is to cut off the spectrum on the halfwidth of gi appr- K 
gives the halfwidht 

Ap = (1 - 2-'/^)po{t) = 0.21poW = /^diffP 
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which is small in comparison with po{g)- So, really, the approximation Qi appr 
can be used here. 

The second way is more close to the iteration procedure from One can 
define Ap by the integral way. One can integrate the approximation Qi appr 
from po up to taking into account that Qi = for p < 0. Then 

Ap = po(t)/4 = AintP 

As far as A^iffP ~ AmtP one can use both these two ways. The integral way 
is more convenient as far as it gives precise asymptotes. 
Now we can suggest approximation 

Gi = fuAppoitf 

and rewrite it as 

Gi=N,{t/A)po{tY 

where Ni{t/4) is the number of the droplets formed until the moment t/4 
and the behavior of G is analyzed in the current moment t. 

When we consider the heterogeneous condensation with essential exhaust- 
ing of centers one can easily note that the function Qi becomes even more 
sharp. So, the previous derivation is suitable here. Certainly, the value of 
Ni(t/4:) has to be calculated with account of exhaustion of the heterogeneous 
centers as it was described in 0, [^. 

Here we can present the last approximation as 

Gi ~ {7]uot-Viit/4:))poitf 
where rji is the number of the free heterogeneous centers of the first sort. 

4 Floating monodisperse approximation 

In PI we were interested in the final parameters of the whole nucleation 
periods and used the monodisperse approximation at t = A2t. It allowed to 
use for A'^i(A2t/4) the following approximation 

iVi(A2t/4) = r]uotil - exp(-5A2t/4)) 

where 

B = fl*nooC*/TT]uot 
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Now we don't want to use the monodisperse approximation only at t = 
A2t. So, we shall act without the last approximation. But now we have to 
use the monodisperse approximation at the arbitrary moment of time t. 

We shall use the variables x, z (see [|l|, 0) and shall investigate the system 
of equations 

= C + Gi + G2 

Gi = fu !\z - x)3exp(ri(C - QIQOi{x)dx 
Jo 

G2 = f2* l\z - xfeM^2{C, - QIQe2{x)dx 
Jo 

9i = exp(- '^^*^°° [ exp(ri(C - C*)/C*)dx) 
Vitot Jo 

02 = exp(- '^^*^°° [ exp(r2(C - C*)/C*)dx) 

V2tot Jo 

where Fj are some parameters (see 0), are the relative numbers of the 
free heterogeneous centers of the given sort. 

In the situation of the strong unsymmetry we can rewrite this system in 
the following manner 

Gi = fu [ {z — xY exp{—TiGi{x)/C^)6i{x)dx 
Jo 

G2 = f2* l\z - xf exp(-r2(G'i + G2)/Qe2{x)dx 
Jo 

9i = exp(- '^^*"'°° / exp(-riG'i(x)/C*)(ix) 

02 = exp(- '^^*^°° / exp(-r2(G'i + G2)/C*)dx) 
V2tot Jo 

The first and the third equations of the previous system form the closed 
system which allows to consider Gi in the second and the forth equations as 
some known value. According to the monodisperse approximation it can be 
presented as 

G,{z) = ^{l-e,{z/4))z' 

where 

ei{z/A) = exp(-E expi-TifuxyAQdx) 
Jo 
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One can simplify the last expression as 



for 
and 

for 

where 
and 



r],{z/A) = r]Mexp{-Ez/A) (3) 

Z/A < Zm 

Viiz/^)=VtoteM-Ei^y/'A) (4) 

J- 1/1* 



Z/A > Zm 

POO 

A = / exp{-x*)dx = 0.905 
Jo 



4C u 



Then the nucleation on the second sort centers can be described by the 
following equations 

G2 = /2* f\z - xf exp(-r2(%(l - 9i{z/4))z' + G2)/Q02{x)dx 
Jo h 



e, = exp{-^^ rexp(-r2(%(l - e^ix/4))x' + G2)/Qdx) 
V2tot Jo h 



V2tot 

We can adopt with a rather high accuracy the following expression for 62 
(the reasons are the same as in the section "Final iterations" in 0) 



9, = exp(-^^ rexp(-r2(4f (1 " 9i{x/4))x' + h.x'/A)/Qdx) 
V2tot Jo E 

or with the help of approximation @, (^) it can be presented in more simple 
form. 

If 2;/4 < Zm then 



9^ = exp{-^^ r exp(-r2(4f (1 - exp{-Ex/4))x' + h.x" lA) IQdx) 
V2tot Jo E 

(5) 
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If 2;/4 > Zm then 

0^ = exp(-:^( T'" exp(-r2(%(l - exp{-Ex/A))x' + h.x'/A)/Qdx + 

r eM-U^a - eM-E{^y/'A))x' + h^x^A) /Qdx)) 

Now we have to calculate the integrals appeared in the last two expres- 
sions. We shall start from the first one. 
Consider (|^). One can see that function 

(j) ^ 12(^(1 - exp{-Bx/A))x' + /2.x74)/C 

is very sharp function. It is more sharp than 

00 = constx^ + constx'^ 

The function (1 — exp(— i?x/4)) is rather smooth in comparison with (j) 
and 00- 

One can note that integrals exp(— a;^)(ix = 0.89 and exp(— a;'^)(ia; = 
0.90 are approximately equal. Both subintegral functions have very sharp 
back front and one can speak about the cut-off in both cases. The approxi- 
mate equality of these integrals means that these cut-off have approximately 
same values. 

We shall define the characteristic parameter Zq by equality 

r2(^(l - eM-EzjA))zl + /2*4/4)/C = 1 (7) 
Then the integral in (|^) can be rewritten as 

rexp(-r9(%(l - exp(-Ex/4))x3 + 
Jo E 

(8) 

f2*x^/A)/Qdx = e{z - Czg)Czg + e{Czg - z)z 

where 

^ roc rOG 

C = -{ exp{—x'^)dx+ / exp{—x^)dx) 
2 JO Jo 
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This representation of the integral transfers into 

= exp(-^^^(e(z - Czg)Cz, + eiCz, - z)z)) (9) 

f]2tot 

Now we shall analyze the integral in (P). The reasons are the same. We 
shall introduce parameter zi by the following relation 

- exp{-E{^f'A))zf + h.zt/A)/C. = 1 
^ J- 1/1* 

Note that we need only one parameter as far as 

1 - exp{-Ex/A) < 1 - exp(-E(-^)i/M) 

for J < Zm and 

1 - exp(-Ea:/4)U./4=.. ^ 1 - exp(-i?(-i^)VM) 

J- ill* 

If zi < Azm then 

/ exp(-r2(^(l - exp(-Ex/4))x3 + /2,x74)/C)rfx > 
Jo h, 

(10) 

r exp(-r2(4f (1 - exp(-E(-i^)VM))x3 + /2*a;V4)/C)rfa; 

"'42m -C/ i l/l* 

and one can analyze only 

= exp(-r2(%(l - exp(-£;x/4))x3 + f2.x^/A)/Qdx 
Jo L 

It was already done in consideration of (^. 
If zi > 4zm then both 

rizm f, 

/i = / exp(-r2(^(l - exp(-Ex/4))x=^ + f2*x^/4:)/Qdx 
Jo h 

and 

l2= r exp(-r2(%(l-exp(-E(-i^)V4A))x3 + /2,xV4)/C)rfa; 
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are essential. 
Then 

h = 4:Zm 

and I2 can be analyzed quite analogously. Namely, we shall introduce Zt from 
equality^ 

U^il - exp{-E{^y/^A))zf + h.zt/A)/C^ = 1 

-c/ i 1/1* 

If Zt is near 4zm then I2 is small in comparison with Ji and there is no 
need to analyze l2- 

If I2 is essential in comparison with Ji one can use the following approx- 
imation: 

l2= {Z- AZm) 

for z < Czt 

h = Czt - iZm 

for z > Czt- 

This completes the approximate analysis of the expression for ^2- 

Some parameters z/, Zm, Zg, zt may coincide but they are conserved in 

order to avoid misunderstanding. 

The main interesting value is ^2(00). The final expressions for this value 

are more simple. They can be directly obtained from the already presented 

ones. 

On the base of 6i{z) one can easily find the number of droplets Ni as 
To find the total number of the droplets one has to put the arguments to 00. 
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